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ABSTRACT 

Globular clusters (GCs) in the Milky Way have characteristic velocity dispersions that are consistent 
with the predictions of Newtonian gravity, and may be at odds with Modified Newtonian Dynamics 
(MOND). We discuss a modified gravity (MOG) theory that successfully predicts galaxy rotation 
curves, galaxy cluster masses and velocity dispersions, lensing, and cosmological observations, yet 
produces predictions consistent with Newtonian theory for smaller systems, such as GCs. MOG 
produces velocity dispersion predictions for GCs that are independent of the distance from the galactic 
center, which may not be the case for MOND. New observations of distant GCs may produce strong 
criteria that can be used to distinguish between competing gravitational theories. 

Subject headings: gravity: theory — dark matter — globular clusters 



1. INTRODUCTION 



Modified Gr a vity ( MOG, iMoffati (f2005l [2006aM : 
iMoffat fc Tothl (|2007bf )) is a fully covariant theory of 
gravity that is based on postulating the existence of a 
massive vector field, cp^. The choice of a massive vector 
field is motivated by our desire to introduce a repulsive 
modification of the law of gravitation at short range. The 
vector field is coupled universally to matter. The theory 
yields a Yukawa-like modification of gravity with three 
constants: in addition to the gravitational constant G, 
we must also consider the coupling constant u> that de- 
termines the coupling strength between the field and 
matter, and a further constant (i that arises as a result 
of considering a vector field of non-zero mass, and con- 
trols the coupling range. In the most general case, these 
constants must be allowed to run. An approximate solu- 
tion of the MOG field equations (| Moffat fe Tothl l2007bh 
allows us to compute the values of fi and u> as functions 
of the source mass. 

MOG was used successfully to describe observational 
phenomena on astrophysical and cosmological scales 
without resorting to dark matter. The theory correctly 
predicts galaxy rotation cu rves (|Brownstein fc M offat 
l2006allMoffat fc Tothll2007bl), the mass and thermal pro- 
files of clusters of gala xies (jBrownstein fc Moffat l F2006b; 
IMoffat fc To thl l2007bD. the merging o f the two clusters 
(Bullet Cluster. lBrownstein fc Moffatj (|2007l )). th e acous- 
tical peaks of the cosmic mic rowave background {Moffat 
I2006bb IMoffat fc Tothl l2007al). th e velocity dispersion of 
satellite galaxies ( |Moffat fc Tothl I2007D . the mass power 
spectrum and the luminos ity-distance relation ship of dis- 
tant Type la supernovae (jMoffat fc Tothl l2007aT l. 

Globular clusters (GCs) in the Milky Way provide a 
particularly interesting case for testing alternate grav - 
ity theories (jBaumgardt et al.lf2005t IScarpa et al.ll2007f L 
such as MO G or Milgrom's Modified New t onia n Dynam- 
ics (MOND. iMfigroml (fl98l : iBekensteinl (l200l L 

GCs at different distances from the galactic center ex- 
perience galactic gravity at varying strengths. The in- 
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ternal gravitational field of a GC also varies depend- 
ing on the mass (typically, 10 4 -10 6 M Q ) and size (typ- 
ically, a few pc to a few ten pc in diameter) of the 
GC in question. Using the characteristic acceleration 
(a = 1.2 x 1CT 10 m/s 2 ) of MOND, for example, we find 
GCs that experience either galactic or internal accelera- 
tion above, or below this value. MOND predicts velocity 
dispersions different from the Newtonian prediction for a 
GC whose stars experience a combined acceleration less 
than ao- 

On the other hand, MOG predicts identical velocity 
dispersions for GCs of the same size and composition, 
regardless of their distance from the galactic center. For 
this reason, GCs provide a unique test by which different 
gravitational theories can be compared. 

2. THEORY 

The bulk properties of GCs, with the possible excep- 
tion of their innermost regions, can be modeled using the 
collis ionless Boltzmann equation (iBinnev fc Tremaind 
Il987| ). from which the statistical properties of the ve- 
locity distribution of stars can be derived. In particular, 
one can derive a formulation for the velocity dispersion 
tensor that, in the isotropic, non-rotating case, reduces 
to a scalar quantity. This quantity can be determined 
using the appropriate Jeans equation. For this calcula- 
tion, one requires knowledge of the distribution function 
(DF) that determines the number of stars in a given re- 
gion of space, and the gravitational potential. We begin 
our discussion with the latter. 

2.1. Modified Gravity 

Our modified gravity (MOG) theory predicts a 
Yukawa-like modification of Newton's law of acceleration 
(|Moffatll2006al : IMoffat fc Totbll2007H ). in the form 



omog 
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(l + a(l-(l + M r)e-^)), (1) 



where Gm is Newton's gravitational cons tant. 

In a ccordance with our recent results (jMoffat fc Tothl 
2007b), the parameters a and [i can now be predicted: 
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where 



v/I7' 



(3) 

G oo ^20G N , (4) 

C[ ~ 25000 M^ 2 , (5) 

C2-6250 MQ /2 kpc _1 . (6) 

For even a large GC, with mass exceeding 10 6 M Q , the 
predicted values are 

a -0.03, (7) 
/j~(160pc) -1 . (8) 

Given the smallness of a and the fact that is much 
larger than the GC radius, it is clear that our theory 
predicts Newtonian behavior for such GCs: 
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For smaller GCs, the predictions are even closer to New- 
tonian values. 

In contrast, the MOND acceleration &mond is given 
by the solution of the non-linear equation 
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where a = 1.2 x 10 8 cm sec 2 . The form of the 
function ji{x) originally proposed by iMilgroml (|1983l ) is 
given by n(x) = x/vl + x 2 ; however, better fits and 
better asy mptotic behavior are ac hieved using fx(x) = 
x/(l + x) (Fa maev fc Binnevll2005D . which yields the ac- 
celeration function 
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Regardless of the form of fi(x) chosen, when the com- 
bined acceleration experienced by stars in a GC is below 
ao, MOND predicts dynamical behavior that is markedly 
different from the Newtonian prediction. 

2.2. The Jeans equation 

In the spherically symmetric, non-rotating case the 
Jeans equation for the velocity disp ersion o~(r) takes the 
follow ing form (see Eq. 4-64a in iBinnev fc Tremaind 

d(yo 2 ) , d$ _ n 



di 
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where r is the radial distance from the GC center, v{r) 
is the number density distribution function, and $(r) is 
the gravitational potential. If v(r) and <i>(r) are known, 
the velocity dispersion can be obtained by direct inte- 
gration. Using a(r) — d^/dr and utilizing the fact that 
lim a 2 (r) = 0, we get 

r — >oo 

oc 

a 2 (r) = i J va{r r ) dr' . (13) 

r 

The observed velocity dispersion of GCs is a function 
not of the actual radial distance r but the projected 



TABLEl 

Properties of GCs studied by IScarpa et al.I (120071) . Data 
for AM 1 and Pal 14 are also included. The distance R g 

FROM THE GALACTIC CENTER, LUMINOSITY Lq IN UNITS OF SOLAR 
L UMINOSITY, AND THE HALF-LIGHT RADIUS Th (PC) ARE SHOWN 
(HARRIS 1996). MASS-TO-LIGHT RATIOS ARE ESTIMATED BY 
FITTING THE VELOCITY DISPERSION USING THE HERNQUIST MODEL, 
EXCEPT FOR AM 1 AND PAL 14, FOR WHICH M/L = 2 WAS FIXED. 



Name 


R g (kpc) 


Lq 


Vh (pc) 


Mq/Lq 


NGC288 


7.4 


3.94 x 10 4 


2.9 


4.38 


NGC5139 


6.4 


1.04 x 10 6 


6.4 


2.79 


NGC6171 


3.3 


5.65 x 10 4 


5.0 


2.54 


NGC6341 


9.6 


1.51 x 10 5 


2.6 


1.50 


NGC7078 


10.4 


3.70 x 10 s 


3.2 


0.85 


NGC7099 


7.1 


7.45 x 10 4 


2.7 


1.51 


AM 1 


123.2 


6.08 x 10 3 


17.7 


2 


Pal 14 


69.0 


6.19 x 10 3 


24.7 


2 



distance R between the GC center and the star being 
observed. The velocity dispersion ctlos(-R) of stars ob- 
served along the line-of-sight (LOS) at projected distance 
R from the GC center is related to a(r) as 



'LOS 



where 
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as can be verified by simple geometric reasoning, 
nating y, we can rewrite (|14|) as 
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2.3. Density Distribution 

Several models exist that can mimic the density dis- 
tribution of a spherically symmetric set of star s . On e 
particularly simple model is that of iHernquistl (|1990h . 
which models the number density of stars as a function 
of radius as 

( \ Nr ° n>7\ 

IHernquistl) = " 1 ; rg, (17) 

27rr(r + ro) 

where N is the total number of stars, and ro is a charac- 
teristic radius. 
Another, similar model is that of iJaffd (|1983f) : 



^Jaffe(^) 



Nrn 



4irr 2 (r + ro) : 



(18) 



Without benefiting from a photometric profile of the 
globular cluster under investigation, there are no a priori 
reasons to prefer one model over another. We found that 
the choice of model does not have a significant impact on 
the conclusions we present; hereinafter, we shall be using 
Hernquist's model consistently, but we note that similar 
results are obtained using alternate number density dis- 
tribution functions. 

3. OBSERVATIONS AND PREDICTIONS 

Velocity di spersion data for sev eral GCs were recently 
published bv lScarpa et al.I (2007). We read velocity dis- 
persion values and their standard devi ations from Fig- 
ures 1-2 and 4 of iScarpa etaTI (|200l . for NGC 288, 
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Fig . 1. — Fitting velocity dispersions obtained from the Jeans equation to globular cluster data (blue error bars from IScarpa et all 
(2007)), using the Hernquist model and MOG or Newtonian gravity (dash-dot green line). 
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Fig. 2. — Predicted velocity dispersion curves for two distant GCs. Dash-dot line (green) is the prediction obtained using MOG or 
Newtonian gravity; dashed (brown) curve is the MOND prediction. In both cases, we used M/L = 2 and we equated the parameter ro of 
the Hernquist model with the half-light radius. 



NGC 5139 (w Centauri), NGC 6171 (M107), NGC 6341 
(M92), NGC 7078 (M15), and NGC 7099 (M30). Some 
of the basic characteristics of these GCs are summarized 
in Table ED 

Using the Hernquist distribution as the number density 
distribution function for a spherically symmetric cluster 
of stars with isotropic velocity dispersion, we obtained 
very good fits to the velocity dispersion data (Figure Q]) . 
These results also yield mass-to-light ratios ranging be- 
tween 0.8 < M/L < 4.4 (Table [TJ, which are typical for 
globular clusters. 

For these results, we used the Newtonian gravitational 



potential. These calculations are consistent with Newto- 
nian theory, MOG (given the smallness of the predicted 
value of the MOG a parameter and the large size of the 
parameter /i^ 1 ), and also MOND, as the GCs in ques- 
tion are located relatively near the galactic center, and 
the galactic acceleration is always greater than oq. 

The possible presence of dark matter does not appre- 
ciable alter these results either. A typical dark matter 
density for the gala ctic halo i s ~ 7 .8 x 10~ 3 M Q /pc 3 
(~0.3 GeV/cm 3 ; see ISumnerl (|2002f) ). a density that is 
much smaller than the globular cluster's stellar mass den- 
sity. 



4 



Moffat and Toth 



The situation is different, however, in the case of 
MOND and globular clusters that are located a long dis- 
tance away from the galactic center. To quote iMilgroml 
(1983): "We are then compelled to conclude that the 
internal dynamics of the open clusters embedded in the 
field of the Galaxy is different from that of a similar but 
isolated cluster." For instance, Pal 14, located at 69 kpc 
from the galactic center, would experience a galactic ac- 
celeration of ~ 2.3 x 10~ n m/s 2 , which is well within 
the MOND regime. As this is a low mass cluster of 
stars, its internal accelerations are also significantly be- 
low MOND's ao, except perhaps in the innermost regions 
of the cluster. 

Two distant clusters, AM 1 and Pal 14, are presently 
the subject of an observational project by Kroupa et al. 
As the absolute luminosity of these GCs is known, using 
a (typical) value of M/L ~ 2 we can obtain a crude 
estimate of their mass, allowing us to apply the Jeans 
equation and derive a velocity dispersion profile using 
both Newtonian and MOND gravity. These predictions 
are shown in Figure [5J 

4. DISCUSSION 

For globular clusters with a mass of a few times 10 6 Mq 
or less, MOG predicts little or no observable deviation 
from Newtonian gravity. This is verified by our demon- 
stration that a simple model, using a spherically sym- 
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metric distribution and no velocity anisotropy, can eas- 
ily reproduce the velocity dispersion profiles of several 
diverse globular clusters with varying mass. 

The predictions of neither Newtonian gravity nor MOG 
depend on the distance from the galactic center. The 
same remains true when dark matter is considered; al- 
though the density of dark matter may be a function 
of distance from the galactic center, at predicted dark 
matter densities, the amount of dark matter contained 
within a GC does not noticeably alter the dynamics of 
the cluster. 

The situation is different for MOND; for a low-mass 
cluster, internal accelerations are below the MOND 
threshold of ao — 1.2 x 10~ 10 m/s 2 , and if the cluster 
is far enough from the galactic center, its galactic accel- 
eration is also below this value. For this reason, distant 
GCs may offer a unique method to distinguish observa- 
tionally between MOG and MOND. 
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